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Abstract 

A number of biological processes, such as embryo development, cancer metastasis or wound healing, rely 
on cells moving in concert. The mechanisms leading to the emergence of coordinated motion remain however 
largely unexplored. Although biomolecular signalling is known to be involved in most occurrences of collective 
migration, the role of physical and mechanical interactions has only been recently investigated. In this paper, 
a versatile framework for cell motility is implemented in-silico in order to study the minimal requirements for 
the coordination of a group of epithelial cells. We find that cell motility and cell-cell mechanical interactions are 
sufficient to generate a broad array of behaviours commonly observed in vitro and in vivo. Cell streaming, sheet 
migration and susceptibility to leader cells are examples of behaviours spontaneously emerging from these simple 
assumptions, which might explain why collective effects are so ubiquitous in nature. This analysis provides also 
new insights into cancer metastasis and cell sorting, suggesting in particular that collective invasion might result 
from an emerging coordination in a system where single cells are mechanically unable to invade. 
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Many biological processes involve concerted cell dis- 
placements across large length scales and time scales. 
These include the most fundamental transformations 
taking place during embryogenesis, such as gastrulation, 
neurulation or vasculogenesis. Common forms of con- 
certed motions are i) cell intercalation, through which 
the tissue undergoes a significant change of shape by cells 
exchanging neighbours (fly germ band extension [I] , ver- 
tebrate gastrulation and neurulation [2p| ), and ii) collec- 
tive migration, during which groups of cells coordinate 
their direction of motion with respect to surrounding 
tissues, travelling reliably large distances in the embryo 
(neural crest cells El, fish lateral line 5 ). Collective 
migration also occurs in adult life, often in association 
with regenerative processes, such as wound healing, or 
pathologies, such as cancer metastasis, which often takes 
the form of groups of cells collectively invading other tis- 
sues 6 . Unravelling the physical and biological princi- 
ples driving collective cell migration is key to understand 
these critical aspects of developments, as well as trigger 
new therapeutic treatments for cancer. The large vari- 
ability in the systems exhibiting patterns of collective 
migration suggests that generic principles are control- 
ling these behaviours 7,8, and calls for an analogy with 



for a directional cue 13,14 . One fundamental assump- 



a larger class of systems 9 . 

More generally, collective behaviours in groups of 
motile individuals have been widely studied over the past 
twenty years, both from a theoretical and experimental 
point of view. Systems considered include animal and 
human populations 10 as well as bacteria 11 or active 
biopolymers 



12 



Flock models have shown that sim- 
ple interactions between the individuals are sufficient for 
coordinating their direction of motion, without a need 



tion in such models is the existence of a local spatial 
coupling, which tends to align the direction of motion of 
neighbour individuals. In the context of cell populations, 
flock models have already demonstrated that a local me- 
chanical coupling is enough to generate streaming pat- 
terns 15,16 and that coordination enhances the sort- 
ing dynamics of heterogeneous populations ujj. Such 
emerging collective effects impose a rethink of the re- 
quirements for large scale cell coordination and in partic- 
ular of how complex cell interactions must be in order to 
ensure robust migration in complex environments. One 
of the limiting factors at this stage is the difficulty to ex- 
perimentally dissect the contributions of mechanical and 
biochemical processes. 



In this paper, a numerical approach is introduced to 
overcome such a difficulty. A range of model studies, 
each highlighting a different aspect of in-vitro or in- 
vivo collective migration, are explored. These include in 
particular the transition from cohesive epithelia to mes- 
enchymal cell populations, the role of population size 
and confinement, the integration of information within 
the population, and the conditions for collective or soli- 
tary cell invasion in a surrounding tissue. Cell interac- 
tions will be purposefully limited to a small number of 
fundamental processes, such as adhesion, incompressibil- 
ity of the cells, and a time scale for the evolution of cell 
polarity. This approach provides, as a result, a generic 
framework from which the role of biomolecular signalling 
can be reinterpreted. 
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1 Results and Discussion 

1.1 Dense endothelial cell monolayers 

The migratory patterns of cultured endothelial cell 
monolayers are well documented, confirming the emer- 
gence of collective behaviours such as streaming and 
large scale velocity correlations 18 - 22 . These experi- 



ments have shown in particular that cell density is one 
of the key parameters controlling collective behaviours. 
In vivo, cells are however usually closely packed, form- 
ing small confined groups within larger tissues [Tp] • The 
onset of group migration then primarily involves varia- 
tions of motile behaviour, adhesive properties or envi- 
ronmental conditions 7 , rather than cell density. The 
first question addressed concerns the role of these physi- 
cal and mechanical parameters in the dynamics of large 
and dense cell populations. The effect of the population 
size and its mechanical environment will be subsequently 
studied. 

The numerical tissue is made of a 2D confluent layer 
of cells on a substrate. Cells can adhere to each other, 
and the resulting cohesion of the tissue is accounted for 
by a membrane tension, J, which is homogeneous across 
the population. J controls here a number of otherwise 
independent properties, such as the adhesion energy, the 
cortical tension, and the amplitudes of membrane fluc- 
tuations 23 . Moreover, all cells have the same volume 
and cannot overlap each other. In this first section, all 
cells also generate the same motile force \i on the sub- 
strate. However, each cell has its own direction of polar- 
ization rii along which this force is produced (|n^| = 1). 
The population lives on a surface with periodic bound- 
ary conditions, such that if a cell leaves on one side, it 
reappears on the opposite. As a result, although the 
system has a finite size, it nevertheless has no boundary 
that could influence cell trajectories. The Cellular Potts 
Model (CPM) [24] , which has already demonstrated its 
broad biological relevance for epithelial tissues, is used 
to implement these rules (see Materials and Methods). 

No specific molecular mechanism for inter-cellular 
communication is considered. Only excluded volume in- 
teractions and constraints of membrane geometry cause 
cells to feel each other or feel any other physical obsta- 
cle. It is however assumed that there is a feedback from 
earlier displacements to the polarization itself, so that 
the cell tries to move towards its direction of polariza- 
tion, but the polarization itself evolves towards the cell's 
net displacement with a response time r (see Materi- 
als and Methods), r acts as a persistence (or memory) 
time for the cell polarization. These assumptions, first 
introduced by Szabo et al. in 2006 in a self-propelled 
model [l5], are consistent with a number of recent obser- 
vations [25 26 , and provide a simple way to model the 
interaction between a cell and its mechanical environ- 
ment. Szabo et al. [16] have used such a model to study 
the dynamics of endothelial cell monolayers. They have 
in particular explored the role played by the persistence 



time in the resulting streaming behaviour and success- 
fully replicated the experimental shape of the velocity 
correlations in homogeneous cell populations using real- 
istic values of the model parameters. This demonstrates 
the general potential of the self-propelled cellular models 
for the study of collective migration patterns. 

The role of the motile force \i on the dynamics of a 
large cell population is first considered, with an empha- 
sis on qualitative transitions, rather than quantitative 
comparisons with experimental data. In the following, 
the persistence time r and noise level are kept constant. 
At low motile forces, no net motion is observed; the pop- 
ulation essentially behaves like a cohesive epithelium, 
with a strong persistence of the neighbourhood relation- 
ships. At higher motile forces, a transition occurs to a 
regime where cells migrate on the substrate and stream 
over large distances (see figure [IJi and Video S2 [35]). 
The spatial correlations of the velocity field quantify the 
structure of the streaming patterns. Considering a cell 
at the origin moving towards the right, the average de- 
gree of alignment between the velocity of this cell and the 
velocity of the cell located at Sr is given by the function 
g(Sr): g = 1 means full alignment, g = — 1 corresponds 
to cell running in opposite directions, and g = to a 
non-correlated situation (see Materials and Methods). 
Figure [TJb) shows 2D color-coded maps of correlation 
functions for three different values of the motile force fi. 

For intermediate values of the motile force (fi « 0.1), 
strong correlations in the velocity field develop, in partic- 
ular along the direction of the cell displacement (stream- 
ing), whereas displacements tend to be anti-correlated in 
the normal direction (meaning that streams of opposite 
direction are running next to each other). Such patterns 
are typically observed in large cell populations T6fl8]|20 
and in-silico |16 . When varying (at constant density) 
the motile force towards smaller or larger values, we find 
however that spatial correlations tend to become more 
localized, and could decay as fast as over a few cell diam- 
eters. The non-monotonic evolution of the correlation 
patterns shows that the dynamics of such a dense cell 
population involves several competing effects that will 
be developed further in the following sections. 

1.2 Emergence of collective motility 

The transition from a non-motile to a motile population 
can be captured by the temporal correlations of the cell 
velocities, or equivalently by quantifying the evolution of 
the mean square displacement of the cells Ar 2 (At, J, /1), 
between two time points t and t + At (see inset of figure 
[2]). For a given pair of parameters (/i, J), the mean- 
square displacement curve is fitted over the time range 
10 < At < 10 4 by the function Ar 2 = DAt^ . The value 
of /3 characterizes the nature of cell displacements over 
that time interval. /3(/i, J) is reported in figure [2] For 
reference, j3 = 1 indicates a stochastic, diffusive like, be- 
haviour. The data show that for each value of J, there is 
a critical motile strength /i c , scaling with J, below which 
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Figure 1: (a) Image of a motile tissue in the steady state, 
overlayed with the corresponding velocity field. Cell col- 
ors are arbitrary, (b) Maps of the velocity correlations 
around a cell migrating towards the right (horizontal axis 
goes from back to front). The unit distance corresponds 
to the cell diameter. 



transport is sub-diffusive (/? < 1), i.e. where cells essen- 
tially never move more than their own size. For /i > /i c , 
migration is however super-diffusive, indicating that the 
direction of motion is persistent over large time-scales, 
several hundred times the memory time (r = 10 in these 
simulations). This regime is consistent with the stream- 
ing patterns discussed in the previous paragraphs. Even- 
tually, as ji increases further, the behaviour progressively 
converges towards a simple diffusive behaviour. At a 
qualitative level, the presence of an optimum for the col- 
lective properties reflects two competing processes taking 
place simultaneously: i) as the motile force increases, the 
coupling between cells also increases, facilitating coordi- 
nation and triggering collective migration. However, ii) 
highly motile cells can also easily disrupt the coordina- 
tion of motile groups by penetrating them. This regime 
of high motile forces will be explored further in the con- 
text of cell invasion (section 1.5). The typical shape 
of the mean-square displacement curves reported in fig- 
ure [2] is in agreement with the theoretical predictions of 
Peruani and Morelli [27] for systems with uncorrelated 
fluctuations for cells speeds and directions of motion; 
speed fluctuations are here related to the random mem- 
brane fluctuations, whereas the polarization evolves over 
larger scales, depending on the cell environment and the 
emerging streaming patterns. These different scales of 
fluctuations are also qualitatively consistent with in vitro 
studies [16], although living systems exhibit more com- 
plex signatures, in particular due to cell divisions |22| . 

The emergence of streaming patterns, observed in 
most self-propelled models, is not in itself surprising. 
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Figure 2: Graph of the diffusion exponent f3 as a function 
of /i/J, for four different values of J (2.5, 5, 7.5 and 10 
- arbitrary units). Inset: average mean square displace- 
ment as a function of time for a range of motility \i and 
adhesion J. Time is expressed in MCS (see Materials 
and Methods) and r = 10 MCS. 



However, interpreted within the context of a cell pop- 
ulation, these results show that i) local coordination can 
arise from purely mechanical interactions, and that ii) 
there is a discrete transition from a static epithelium to 
a streaming population. This behaviour is reminiscent of 
the Epithelial-Mesenchymal Transition (EMT) in living 
systems, in which cell behaviour changes rapidly from 
non-migratory to migratory. The causes of such a tran- 
sition in vivo are still highly debated 28 29 . Although 



a large number of oncogenes are involved, the EMT does 
not seem to be associated with a proper lineage switch. 
Figure [2] shows that a number of continuous physical 
quantities, such as the motile force or the cohesion of 
the cellular layer, control a rapid transition between two 
discrete dynamical states, from a static tissue to a highly 
mobile population of cells, and vice- versa. It appears 
that, even in such a simplistic model, a significant change 
in a tissue phenotype can originate from minute modula- 
tions of either the cells' behaviour or their environment, 
without a need for a cell-lineage switch; an extrapolation 
to more realistic 3D situations might explain the practi- 
cal difficulty in finding robust determinants for the EMT, 
in particular based on gene expression. 

1.3 Emerging consensus and size effects 

In vivo cases of collective migration usually involve small 
groups of cells with highly coordinated directions of mi- 
gration, for instance during the collective invasion of 
carcinoma cells, or in a number of developmental pro- 
cesses such a neural crest cell migration or lateral line 
primordium (migrating epithelium). In order to probe 
further the effect of the population size, systems contain- 
ing from 9 to 1600 cells with periodic boundary condi- 



3 



tions have been studied. The substrate size was changed 
accordingly so that the cell density remained constant 
(see Video S3 [35]) . The degree of global coordination 
is quantified as the mean velocity across the whole pop- 
ulation normalized by the mean cell speed ((v) / (\v\)). 
This corresponds to an order parameter, taking values 
from (no order) to 1 (full coordination or sheet migra- 
tion). Figure |3^i shows how this parameter depends on 
the motility, for various values of J. For each set of pa- 
rameters, one can identify a typical length scale A p (/i, J) 
corresponding to the largest system size at which global 
coordination spontaneously arises. As expected from 
the previous results, the ability to align is maximum 
when the system is at the onset of collective migration 
(/i > /i c ), i.e. when temporal and spatial correlations 
are at their maximum. Sheet or group migration ap- 
pears therefore as an intrinsic behaviour of groups of 
motile cells: for a small cell population, an increase of 
cell motility leads to a spontaneous transition from a 
non- motile epithelium to a migrating sheet. These re- 
sults might shed light on recent observations of sheet 
migrations in vivo, for instance in the fish lateral line 
primordium pi. The model predicts in particular that 
a progressive increase of the traction force (or decrease 
of cohesion energy J) in an epithelium leads first to a 
regime of high coordination, where sheet migration is 
relatively straightforward to achieve. 

1.4 Leader cells and integration of exter- 
nal cues 

Coordination alone is however not enough to prescribe 
any particular direction of motion for the population. 
Setting a directional preference for a group of cells re- 
quires some level of interaction with the environment, 
for instance to sense a gradient. In the context of collec- 
tive migration and invasion, only a small proportion of 
competent (leader) cells in the population seems to be 



required in order to induce directed motion 19 , 20 , 32 



However, how these leaders might steer the whole pop- 
ulation remains unclear. In this section, physical mech- 
anisms by which leader cells could influence large scale 
migratory patterns are explored. 

Within our framework, leader cells would be cells 
whose polarity is set by an external cue, rather than 
through a feedback from their previous displacement; 
their direction of polarity n is therefore kept constant, 
aligned along a fixed arbitrary direction, identical for all 
leader cells. These cells are otherwise indistinguishable 
from the other cells (including same motile force \i and 
tension J). By contrast with live situations where leader 
cells tend to be located at free boundaries, a small pop- 
ulation of such leader cells is here scattered uniformly 
in the tissue (see figure [3^ and Video S4 [35]). This 
choice allows us to study specifically how information 
propagates through the population; modelling a free sur- 
face would introduce many additional parameters which 
would make the process less tractable at a qualitative 
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Figure 3: (a) and (b): Heat maps of the order parame- 
ter (v) I (\v\) as a function of (a) /i and the system size 
d g (number of cells: N « d 2 g ) or (b) \i and the typi- 
cal distance between leaders d\ (in cell diameters). Data 
for three different values of J is presented, (c) Exam- 
ple of a tissue with a few leader cells (with pink/orange 
tone) whose polarity is constant and directed towards 
the right, (d) Sketch of the curve X g (fi) and its qualita- 
tive relationship with the different regimes of migration. 
For a given length scale d associated with a constraint 
(distance between leaders, distance between boundaries, 
or number N of cells in the group (d = y/N)), three 
regimes can be defined as \i increases: epithelium, sheet 
migration, or uncoordinated. 



level. 

Considering large populations (too large to spon- 
taneously coordinate), the coordination parameter 
(v) / (\v\) is now calculated as a function of the cell 
motile force \i and the distance d\ between leader cells 
(see figure ^jp) . This distance is now the most relevant 
length scale, related with the leader density pi « d^ 2 . 
The larger the leader density (or the smaller the dis- 
tance between leaders), the stronger is the resulting co- 
ordination. However, as expected from systems of self- 
propelled particles [30], the susceptibility of the popu- 
lation to leader cells strongly depends on the collective 
aspect of the dynamics; when /i « /i c , 1% of leader cells 



4 



are sufficient to significantly coordinate the whole cell 
population. As previously, one can define a length-scale 
for each set of parameters, A/(/i, J), corresponding to 
the largest distance between leaders that can induce a 
global motion of the population. A comparison between 
the figures [3^i and ^jp shows that the susceptibility to 
leader cells strongly correlates with the intrinsic ability 
of the population to coordinate their direction of mo- 
tion. A population will respond well to leader cells if 
d\ < A/(/i,J) « \g(jjL,J). The population's response 
to leader cells is therefore primarily controlled by its 
own internal dynamics, encompassed by the length-scale 
\ g (/j,,J). Any constraint that acts over a length scale 
d smaller than A^(/i, J) is susceptible to trigger global 
coordination, as summarized in figure [3]i. In particular, 
any persistent directional bias (leader cells, anisotropic 
boundaries, etc) could steer the overall migration, as 
long as coordination is already ensured by the collective 
dynamics of the group. 

Another implication of this model is that neither 
specific mechanisms of communication between leaders 
and normal cells, nor a long range mechanical coupling 
through the substrate 21 , are required for large scale co- 
ordination; the same short-range mechanical interactions 
are enough to serve that purpose. The idea that sim- 
ple physical/mechanical interactions could indeed play 
a role in the leader /follower relationship has been men- 
tioned [5{|9|[l9], and the results presented above demon- 
strate that such a scenario is indeed plausible. The be- 
haviour observed is moreover in very good agreement 
with the in vitro studies of Vitorino and Meyer 20 who 



have shown that i) local coordination within 2D popula- 
tions of endothelial cells is independent of the presence 
of directed migration in the population, and that ii) only 
a very small proportion of leader or "pioneer" cells are 
required to steer the whole population of cells. 

1.5 Tissue invasion 

Collective forms of tissue invasion not only involve the 
coordination of a population of cells, but also its ability 
to penetrate a surrounding matrix or tissue. A number 
of physical and biological factors are known to influence 
collective invasion 6 , but a general picture remains chal- 
lenging to extract at this stage. For instance, both the 
intrinsic behaviour and the environment of cancer cells 
seem to determine the collective or individual nature of 
the invasion process. In this part, the framework in- 
troduced above is adapted to study cell invasion in a 
passive surrounding tissue. The tumour cells are essen- 
tially a population of motile cells enclosed by a tissue. 
The latter is modelled for simplicity as an epithelium 
of non-motile cells with the same mechanical properties 
(encompassed by J) as the tumour cells, so that adhe- 
sion based sorting effects will not occur [24]. J is kept 
constant in what follows; its role would simply be to shift 
the transitions as observed on figures [2] and [3] The effi- 
ciency of invasion as a function of the motile strength of 
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Figure 4: (a) Diffusion exponent (3 of a single motile cell 
in a tissue a non-motile cells. Left inset: Mean square 
displacement of single motile cells. Right inset: sample 
image of a single motile cell (blue) in the tissue, (b) 
Graph of the total invasion distance A as a function of 
distance travelled by the motile cells (scaled time), for 
various motile forces (J = 5). 



the tumour cells is now characterized. 

The minimal motile force ji s required for the migration 
of a single cell in the passive tissue is first determined (see 
Video S5 to S7 [35]). Figure[4^ shows the evolution of the 
mean square displacement of individual motile cells in a 
non-motile tissue. ji s « 0.225 corresponds to the cross- 
over between a sub-diffusive behaviour (cell trapped in 
the tissue) and a migratory behaviour. This threshold 
is here larger than the motility \i c required for a motile 
population to undergo an EMT. 

A tumour is now considered, modelled as a disc of 400 
tumour cells (motile force \i > 0) surrounded by thou- 
sands of non-motile tissue cells. The invasive behaviour 
of the tumour is quantified by a scalar quantity A(/i, t) 
corresponding to the sum, for all the cells which have 
left the tumour, of their radial distance from the initial 
tumour boundary; it captures both the distance trav- 
elled and the number of cells migrating away. The time 
evolution of this quantity for a range of motile strengths 
is reported in figure (4Jd. For low cancer cell motility 
(/i < /i c ), cells are not migrating and no invasion oc- 
curs. For high motility (/i > /i s ), cells are motile enough 
to escape individually, leading to a very efficient inva- 
sion. Cells eventually fully disperse in the tissue, and 



5 



A reaches then a plateau due to the finite size of the 
system. Although individual cells cannot invade when 
/i < a significant amount of invasion however takes 
place in the range fi c < fi < In this regime, groups 
of cells, forming transient streams of 5 to 10 cells, are 
frequently seen leaving the tumour together, forming 5 
to 10 cell long trains (see inset in figure and Video 
S8 & S9 citeSupMat). This behaviour is explained by 
the fact that coordinated motile cells generate together 
larger forces and can penetrate the surrounding tissue; 
it is therefore analogous to the "tug of war" process ob- 
served in vitro by Trepat et al. [3l]. Collective invasion 
is favoured in this regime because single cell invasion is 
mechanically impossible, but small groups of cells can 
nevertheless coordinate and join forces. 

Although real tumour invasions are far more complex 
than the system studied here, the typical patterns of col- 
lective invasion only require straight-forward conditions: 
(i) an ability to coordinate cell migration (through the 
mechanisms suggested above, or otherwise) and (ii) a 
mechanical environment (surrounding tissue, or extra- 
cellular matrix) that can resist single cell invasion but 
fails if the force is slightly higher. This might explain 
why a particular cell type could exhibit different invasion 
strategies depending on the properties of its direct envi- 
ronment (represented here by /i s ), as often observed [6]. 
This also suggests that having a specialized cell at the tip 



of the invading group 6 32 is not necessary for coordina- 
tion, but might be mechanically critical. In the common 
case of fibroblast led invasion, collective invasion might 
for instance proceed by i) a track being created in the 
matrix/surrounding tissue by the competent fibroblast 
cells [32] and ii) cancer cells streaming in the new space, 
strongly coordinated due to confinement in the track. 

1.6 Dynamic sorting 

The populations of cells considered in this paper all have 
the same membrane tension J, in order to avoid addi- 
tional surface tension effects leading to cell sorting. In- 
deed, depending on the relative affinity of the different 
cell populations, mixing might be enhanced or prevented, 
interfering with the qualitative behaviour studied in this 
paper. In fact, invasion and sorting are strongly related 
phenomena: a poor ability to invade is analogous to a 
strong tendency to segregate. 

To explore further this analogy, cell sorting is now 
studied in a system containing an equal proportion of 
motile and non-motile cells, initially arranged in a chess- 
board pattern. Cell sorting is monitored by calculating 
r(t, /i) = timm I^mn "j where umm is the number of con- 
tacts between motile cells, and unm is the number of 
contacts between motile and non- motile cells. A large 
value of T indicates a strong sorting effect, whereas a 
random distribution would lead to T = 1. The figure 
[5|a) shows the evolution of T as a function of time, for 
different values of the motile force ji. For ji < /i c , no sort- 
ing occurs since cells are not moving (static epithelium 



state). However, for motile forces just above /i c , sorting 
occurs, eventually leading to the formation of massive 
clusters of non-motile cells, surrounded by large streams 
of motile cells (figure |5|b) and video in supplementary 
material [35]). As the motile force further increases, sort- 
ing becomes nevertheless less and less efficient, as motile 
cells are becoming invasive and easily migrate through 
the groups of non- motile cells. 

Earlier theoretical work on cell sorting has focused on 
adhesion differences to explain the separation of cell pop- 
ulations. Although cell motility is known to influence the 
dynamics of adhesion based segregation [17] , its possible 
role as a driving force remained unexplored. The results 
presented here demonstrate that differences in motility 
can be sufficient to drive the sorting of cell populations, 
even without any significant difference in adhesive prop- 
erties. A distinctive signature of such a process is the 
formation of streams of motile cells, and islands of non- 
motile groups. This prediction remains to be supported 
by experimental data. 
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Figure 5: (a) Evolution of T(t) (the number of contacts 
between two motile cells divided by the number of con- 
tacts between motile and non motile cells) for different 
motile forces /i ( J = 5). (b) Segregating tissue at t = 10 6 
MCS with fi = 0.125 and J = 5. 
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1.7 Discussion 

1.7.1 A unified picture for collective migration 

The practical importance of the ideas introduced here 
relies most of all on their ability to provide, with par- 
simony, a unified picture of the patterns of collective 
migration. This work demonstrates that a non specific, 
mechanical coupling and a persistence time in the po- 
larization direction of the cells are enough to reproduce 
a wide array of cellular behaviours commonly observed 
in vitro and in vivo. It reproduces the main features of 
the in vitro behaviour of endothelial cells away from a 
wound, and explains the mechanisms of large scale coor- 
dination when a wound is present. It shows that sheet 
migration is a robust feature of small groups of migrating 
cells, and predicts that the length scale required for sheet 
migration is related with the dynamics of the population 
in an unconfined condition. Finally, it provides guidance 
to interpret cell invasion processes and how these depend 
on matrix properties. 



malignant fibroblasts 34 . One might question to begin 
with why normal fibroblasts would need contact inhibi- 
tion of locomotion, and a possible answer might lie in the 
phenomena of dynamic segregation. As demonstrated by 
figure |5j unless specific mechanisms are in place, motile 
fibroblasts would have a spontaneous tendency to stream 
and segregate rather than disperse in the tissue. CIL 
might therefore be a necessary behaviour to avoid segre- 
gation. This could also explain why changes in CIL can 
influence the emergence of collective effects in cell pop- 
ulations HI, and in particular help tumour cell invasion. 
If this model suggests generic mechanisms for collective 
behaviours, it also shifts the biological questions toward 
the understanding of how populations of active cells can 
accurately control their dynamic state in critical parts 
of development and tissue homeostasis. A number of 
current signalling mechanisms might have emerged from 
such evolutionary constraints. 

2 Materials and methods 



1.7.2 Dynamic transitions 

The dynamics of a cell population has been primarily 
presented as a function of the motile force, for the sake 
of simplicity. However, the critical values \i c and \i s at 
which the qualitative transitions occur do depend on the 
other physical and biological parameters involved in this 
model. For instance, /i c increases with J (see figure |5| 
and decreases with the persistence time r. As a result, 
transitions can be triggered by a variation of adhesion, 
cortical tension or response time as much as by a change 
in motile force fi. The concept of dynamical transition, 
though common in physical systems, has a number of 
additional implications in a biological context. The re- 
sults presented here show that minute variations of cellu- 
lar and environmental parameters can trigger transitions 
between well defined phenotypes at the population scale 
(epithelium/sheet migration/mesenchyme, coordinated 
motion/disordered motion, no/collective/individual in- 
vasion). This might explain why such changes in pheno- 
type are often difficult to associate with a proper lineage 
transition (or a change in gene expression) and how they 
depend on external parameters. 

1.7.3 Biomolecular signaling 

Highly organized behaviours are emerging without the 
need for specific cues of inter-cellular communications; 
this raises a number of questions concerning the role of 
known signalling pathways identified in association with 
individual or collective cell migration 33 . Cell-cell sig- 



nalling might be seen either as a complementary (or re- 
inforcing) mechanism, or as a way to compensate for col- 
lective effects. While any conclusions at this stage would 
be speculative, one might already comment on a few well 
established facts. Malignant fibroblasts have a reduced 
contact inhibition of locomotion (CIL) compared to non- 



2.1 Numerical model 

These simulations are implemented using the Cellular 
Potts Model (CPM) 24 . The CPM is a lattice model: 
the state of the tissue is discrete, represented by an image 
of the system where the value a(k) of each pixel location 
k codes for the identity of the cell covering that location. 
Cells can a priori adopt any shape on the lattice. The dy- 
namics is introduced by minimizing an energy function 
of this state and other physical and biological param- 
eters of the system. Adhesion and cortical tension are 
accounted for by a single parameter J that provides an 
energy cost per unit of membrane length (in 2D) between 
cells. The cell volume is constrained to a reference value 
^o, with a compressibility The energy function is 

minimized by randomly choosing one pixel at a time and 
testing if the energy can be lowered by transferring that 
pixel to a neighbouring cell. Fluctuations are introduced 
as part of this process, with a typical scale represented 
by a parameter T. Time is expressed in Monte Carlo 
Steps (MCS), where 1 MCS corresponds to an average 
of one iteration per pixel over the whole lattice. 

2.2 Motile force 

Each cell generates a motile force along its polarization 
direction n^, with a strength \±i. These forces are here 
forces between each cell and its substrate. Interactions 
between cells being already accounted for by the mem- 
brane tension term. The energy function used in the 
CPM accounts for motility by adding, for each cell z, a 
time dependent energy gradient along n^. These essen- 
tially act as sources of energy that can drive cell motility. 
The function can be written: 

E = Y>i-K{vi -v ) 2 + ^k,iJ • (l - 5a(k),a(i)) - ^ifiiTii-ri 

(1) 
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where sums over i are sums over all cells, and pairs (k,l) 
represent neighbouring pixels. S a ^\ a m is 1 when the 
two pixels belong to the same cell, and otherwise. The 
first term enforces the volume constraint, the second the 
membrane tension, and the third drives the motility. 

Dissipation in this model occurs due to the energy fluc- 
tuations, resulting in an effective viscous friction force 
between the cells and the substrate. The mean cell speed 
within the epithelium is mostly a linear function of \±j J 
for the range of motile forces used here. This is due to 
the fact that J controls the amplitude of the membrane 
fluctuations, which in turn controls the dissipation be- 
tween the cell and the substrate (fluctuation-dissipation 
theorem). In such a Potts model, the dynamics is there- 
fore implicitly overdamped. 

2.3 Dynamics of the cell polarity 

The direction n of the motile force is determined by a 
feedback from its earlier displacements (unless specified 
otherwise, for instance concerning leader cells), n is ori- 
ented along the mean velocity of the cell over its past r 
time-steps: 

( V i)[t-r,t] 



rti(t) 



(Vi) 



[t-r,t] 



(2) 



The parameter r characterizes the time-scale at which 
cell polarity evolves and integrates changes in the me- 
chanical environment, such as contact with boundaries 
or with other cells. 

2.4 Parameters used 

The simulations presented here have been performed 
using next-nearest neighbours interactions in the en- 
ergy calculation, a preferred cell area (^o) of 400 pix- 
els (20x20), an energy fluctuation scale T=2.5, a com- 
pressibility ft -1 =0.5 and a memory time r=10 MCS. The 
results of figure [3] have also been reproduced using next- 
next nearest neighbour interactions, using up to 6400 
cells. The initial state of the system corresponds to a 
regular tessellation of the substrate, with each cell start- 
ing with a random orientation of its polarization. 

2.5 Spatial correlations of the velocity 
field 

The analytical expression of g(Sr) is: 

(v (r) • v (r + R (r) Sr)) r 



g(Sr) 



(v 2 ) 7 



(3) 



where R(r) is the transformation rotating v (r) along 
the horizontal x axis. 
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